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Abstract 



We report on a study of a bilayer two-dimensional electron gas at Landau 

level filling factor u = 2. The system exhibits both magnetic and spontaneous 

interlayer phase coherence broken symmetries. We propose a 3-parameter 

Slater determinant variational wavefunction which describes the ground state 

over the full range of bias potential (Ay) Zeeman coupling (A z ) and interlayer 

tunneling amplitude (At) strengths. Broken symmetry states occur inside a 

volume in this three-dimensional phase diagram near the A 2 = Ay + A 2 

surface. We have obtained analytic results for the intersections of the phase 

diagram with the At = 0, A z = 0, and Ay = 0, planes and show that the 

differential capacitance of the bilayer system is singular at the phase boundary. 
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I. INTRODUCTION 

The interaction physics of two-dimensional electrons in the quantum Hall regimeS is 
enriched by the macroscopic degeneracy of zero-width Landau bands. Among the conse- 
quences of electron-electron interactions in the quantum Hall regime are broken symmetry 
ground states that are especially robust at integer Landau level filling factors. In a single- 
layer two-dimensional electron gas, the ground state at Landau level filling factor v — 1 is 
a strong ferromagnetH with total spin quantum number S = N/2 where N is the number 
of electrons. When mixing between different orbital Landau levels is neglected, the exact 
ground state is known and consists of a fully occupied set of lowest-Landau-level (LLL) 
orbitals sharing a common spinor characterized by an arbitrary! spin orientation. This bro- 
ken symmetry ground state is conveniently described using the lexicon of the Hartree-Fock 
mean-field approximation, which happens to be exact in this instance. In this language, 
the broken symmetry is a consequence of exchange interactions which spontaneously spin- 
split the orbital Landau levels. The Stoner criterion for ferromagnetism in itinerant-electron 
mean-field theory, Ig{tF) > 1 where I is the exchange integral and g{eF) is the density-of- 
states at the Fermi energy, is always satisfied because the density-of-states is infinite for a 
zero-width band. 

For bilayer quantum Hall systems, described schematically in Fig. |I], the Hartree-Fock ap- 
proximation ground state at v = 1 breaks both spin-rotational invariance and the individual- 
layer charge conservation symmetryil, by fully occupying an orbital with an arbitrary! spin- 
orientation and spontaneously developed phase coherence between top and bottom layers 
even in the absence of interlayer tunneling and Zeeman coupling. Bilayer systems are richer 
still for Landau level filling factor v — 2. In the Hartree-Fock approximation, two orbital 
Landau levels are fully occupied, the spin and bilayer state of each described by a four- 
component spinor. In realistic circumstances electron-electron interactions can be more 
important than the one-particle interaction terms corresponding to Zeeman coupling, bias 
potential, or inter-layer tunneling energy in determining the form of these spinors, opening 
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the possibility for broken symmetries in the ground state. In previous workO it has been 
established that the Hartree-Fock ground state at zero bias potential and weak Zeeman 
couplings is a canted antiferromagnet, with spin ordered moments having opposing tilts 
away from the Zeeman field direction in opposite quantum wells. Theoretical predictions 
for the canted antiferromagnet state are in good agreement with the findings of inelastic 
light scattering experimentsil. Although, unlike the single layer case, the Hartree-Fock ap- 
proximation does not yield the exact ground state in bilayers, it can be improved upon 
systematicallyB and should reliably indicate the types of broken symmetry phases which 
occur. In this paper, we use the Hartree-Fock approximation to describe the dependence of 
the ground state of v = 2 bilayer quantum Hall system on a bias voltage which moves the 
two layers away from density balance. A bias potential is readily applied in situ and, as we 
show, can be used to tune the system either in to or out of the broken symmetry region 
of the phase diagram. We point out that the differential capacitance of the bilayer system, 
which is relatively easy to measure experimentally, will show singular behavior at the phase 
boundaries. In Section II we introduce a four- component spinor representation of bilayer 
quantum Hall ferromagnets at v — 2. We propose a general 3-parameter Slater determinant 
variational wavefunction which accurately describes the ground state over the entire phase 
space. We present an explicit expression for the ground state energy functional and the four- 
dimensional Hartree-Fock Hamiltonian whose eigenvalues can be used to estimate the gap 
for charged excitations. Sections III, IV and V report analytical results for the cases of zero 
Zeeman coupling, zero bias potential, and zero interlayer tunneling, respectively. Section III 
and V also include analyses of the singularity which occurs in the bias-potential-dependent 
internal capacitance at the order-disorder phase transition. In Section VI we report numeri- 
cal results for the dependence of several quantities which are open to experimental study on 
model parameters. Finally in Section VII, we conclude with a discussion of the dimension 
of the order parameter space. 
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II. HARTREE-FOCK VARIATIONAL WAVEFUNCTIONS 



Bilayer models for quantum well electron systems apply whenever the growth-direction 
degree of freedom of the electron is adequately described by specifying the well in which an 
electron resides, or equivalently whenever the two lowest electric subbands are well separated 
from higher subbands. We adopt the bilayer model throughout this article. The single- 
particle Hilbert space is then a direct product of the two-dimensional orbital continuum 
space, the two-levels associated with the electronic spin degree of freedom, and a second two- 
level system which constitutes the growth-direction layer degree of freedom. It is convenient 
and has become a standard notation to describe orbital two-level systems using a pseudospin 
language. In this paper we will treat the layer degree of freedom as a pseudospin and use 
a convention in which an electron in the top layer is said to have pseudospin up while an 
electron in the bottom layer is said to have pseudospin down. The magnetic field is assumed 
to be sufficiently strong that only 2D electron orbitals in the LLL can be occupied. We 
further assume that translational invariance is not broken in the ground state. We seek 
the single-Slater-determinant many-particle wavefunction which minimizes the expectation 
value of the Hamiltonian. It follows from the above that it has the form 



Here c\ x is a creation operator for a Landau gauge LLL orbital. The index k specifies the 
spin and pseudospin state; in our convention k — 1 for an up-spin electron in the top layer, 
k = 2 for a down-spin electron in the top layer, k = 3 for an up-spin electron in the bottom 
layer and k = 4 for a down-spin electron in the bottom layer. Translational invariance 
requires that the coefficients z l k are independent of the label X which specifies the guiding 
center state within LLL. At v — 2, the index i = 1,2 and the ground state consists of two 
filled Landau levels specified by four-component spinors z 1 and z 2 , respectively. This is an 
enormous simplification but, at first sight, still leaves us with a sixteen-dimensional space 
specified by eight complex numbers in which to search for a minimum. The considerations 




(1) 



i,X fc=l,4 
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we use to reduce this dimension are dependent in part on the form of the Hamiltonian which 
we now specify. 

The model we examine allows for external bias potential, tunneling between the quantum 
wells, and Zeeman coupling, in addition to electron-electron interactions. We write the 
single-particle Hamiltonian in the form 

h° = -(A v /2)r 2 - (A t /2)r x - (A z /2)a x . (2) 

We have introduced the notation r a for the Pauli spin matrices acting on the pseudospin 
degree of freedom and a a for the those acting on the spin degree of freedom. In the current 
context these are four-dimensional matrices whose explicit form is easily constructed from 
our convention for the k indices specified above. The parameter Ay represents the electric 
potential drop between the two quantum wells, due to electric field created by charges 
external to the bilayer system, A t /2 is the amplitude for an electron tunneling between 
top and bottom quantum wells which can! be measured by studying weak-field Schubnikov- 
deHaas oscillations, and A z is the single-particle spin-splitting due to Zeeman coupling. The 
single-particle terms yield four macroscopically degenerate Landau levels whose energies are 



E° - ±- 

-^1,2,3,4 ~ ^2 



A? + A 2 V ±A Z 



(3) 



Broken symmetry states for v = 2 occur only when the relative and absolute strengths of the 
three single-particle terms is such that the gap between the two lower energy eigenvalues 
and the two higher energy eigenvalues is smaller than or comparable to electron-electron 
interaction strengths. 

The many-particle Hamiltonian whose expectation value we minimize is 

H= J] C l',xK',k C k,X + -X Yl Y c kl,Xl> c k2,X2> c kl,XlCk2 y K2[(X 1 ,X 2 \V + \Xi,X 2 ) 

k',k,X kl,k2 X 1 ,X 2 ,X' 1 ,X' 2 

+ r^ kl r^ k2 (X' 1 ,X' 2 \V_\X 1 ,X 2 )}. (4) 

Here V± = (Vs ± Vb)/2 and Vs and Vd are the 2D interactions between electrons in the 
same layer and in different layers. The qualitative physics of bilayer quantum Hall systems at 



v = 2 is controlled by the property that Vd is weaker than Vs. Note that this model explicitly 
drops the, normally small, terms in which interactions scatter electrons between layers. The 
model also assumes that the two quantum wells in the bilayer system are identical, usually 
true of in experimental systems. This assumption can easily be relaxed and does not change 
any qualitative physics. 

The properties of orbital statesi within a Landau level can be used to derive a general 
expression for the variational energy of the Slater determinant trial wavefunctions: 

1 2 

e = — l A v (z l \r z \z l ) + A t (z l \T x \z l ) + A z (z l \a x \z l ) 

2 i=i 

2 

+ (H{z i \r z \z i ){z j \T z \z j ) - F^{z i \T z \z j ){z j \r z \z i )) ] (5) 

Here e is energy per Landau level orbital. (The number of orbitals in a Landau level is 
= AB/Qq where A is the cross-sectional area and $o is the electron magnetic flux 
quantum.) In Eq.[ |5| we have adopted Dirac notation for the two four- component spinors 
and introduced the three interaction energy parameters which appear in the Hartree-Fock 
theory of bilayer quantum Hall systems: 

H=(2nf)- 1 V-(q = 0), (6) 

and 

F± = J ^- 2 exp(-qH 2 /2)V ± (q). (7) 

The parameter H characterizes the electrostatic energy associated with charge transfer be- 
tween the wells, while F + and _F_ are exchange energies associated with the sum and differ- 
ence of intralayer and interlayer interactions. Note that F + contributes only a constant to the 
energy so that the ground state depends only on the three single-particle parameters and on 
H and _F_. In the limit of infinitely narrow quantum wells and for the LLL, H = (e 2 /e£)d/2£ 
where e is the dielectric constant and £ = (hc/eB) 1 ^ 2 is the magnetic length. The interaction 
parameters depend in general on details of the bilayer system geometry and on the index 
of the orbital Landau level at the Fermi energy. The values of these parameters are easily 

6 



evaluated once the geometry has been specified. We assume in what follows^ that H > F_, 
an inequality which is always satisfied when the interacting 2D electron layers are identical. 

Extrema of the Hartree-Fock variational energy functional satisfy the Hartree-Fock single 
particle equations. For the present problem, the spinors we seek are the two lowest energy 
eigenstates of the four-dimensional matrix 

h HF = h o_ F+p + Hr z [Tr(jrr z )] - F_t z P t z (8) 

where p is the four-dimensional density matrix which must be constructed self-consistently: 

Pkl , k = jZ{b k 'W){^\b k ) (9) 

8=1 

Here b k denote the four four-component basis spinors (6 1 = (1, 0,0,0), 6 2 =(0, 1,0,0), 
6 3 =(0, 0,1,0), 6 4 =(0, 0,0,1)) corresponding to label-/c spin and pseudospin states as introduced 
in Eg. ([!]). The Hartree-Fock variational state can be determined by finding solutions of the 
Hartree-Fock equations or directly by minimizing the energy functional fl5|) with parameter- 
ized spinors. The former approach is usually more convenient when numerical solutions are 
required and has the advantage of yielding an estimate of the quasiparticle excitation gap, 
which is important in interpreting transport experiments^. The latter approach is often 
more convenient when the calculation can be done analytically. In this paper we will try to 
proceed as far as possible analytically before resorting to numerical calculations to illustrate 
and extend the analytical results. 

We start from a variational ansatz for the two four-component spinors z 1 and z 2 ap- 
pearing in the v = 2 translationally invariant trial wavefunction (Eq.[ 0). (A more complete 
discussion of the most general translationally invariant v = 2 wavefunction and of the num- 
ber of parameters it can carry will be given in the concluding Section VII.) We have checked 
numerically to ensure that our ansatz, characterized by three variational parameters, is 
sufficiently general to capture the minimum energy state for any values of the parameters 
which appears in the one-body Hamiltonian (0), namely, the bias potential, the Zeeman 
coupling and the tunneling amplitude. Our wavefunctions are specified by the following two 
four- component spinors: 



z 1 = (Cc^Ss^Cs^-Sci) 

z 2 = (-Ss h C Cl ,Sc h C Sl ) (10) 

Here C = cos(/9/2), S = sin(0/2), c a = cos(x<r/2), and s a = sin(x cr /2) with a = t ? I- The 
variational energy of these wavefunctions is 

e = —Ay 008(6*) cos(x+) cos(x-) — \ cos(9) sin(x+) cos(%_) — A z sin(#) sin(x_) 
-F+- F_[sin 2 (#) sin 2 (x-) + cos 2 (#)(cos 2 (x + ) cos 2 (x-) + sin 2 (x+) sin 2 (x-))] 
+ 2H cos 2 (fl) cos 2 (x+) cos 2 (x-) (11) 

Here x± = {X{ =t Xt)/2- We have not been able to minimize this energy expression analyti- 
cally for the case when all three parameters Ay, A t and A 2 are non-zero. In the following 
sections we present analytic results for the cases in which one of the three parameters is set 
to zero. 

The density matrix corresponding to our variational wavefunction is 

2p = 1 + cos(6 ) ) cos(x + ) cos(x_)r z + sin(6 ) ) sin(x_)cr x ' 

+ cos(#) sin(x+) cos(x-)r x + sin(0)(cos(x_) - cos(x+)K't x 

- cos(6>) cos(x+) sm(x )o z T x + cos(6>) sin(x+) sin(x-)cr z r 2 (12) 

Expectation values of spin and pseudospin operators and their products can be read off 
from Eq. ([12]) by using the familiar properties of Pauli matrices. For example, the difference 
between the z direction projection of the spin in the top and bottom layers, (r z a z ) = 
Tr(pr z a z ) which was identified as the order parameter in the work of Zheng et a/.§, is given 
by 

O zz = 2 cos(0) sin(x+) sin(x_). (13) 

This order parameter is like that of a canted antiferromagnet; indeed this is the rubric used 
by Zheng et al. to describe the broken symmetry state. We will show that 

O xz = -(rV 2 ) = 2 cos(#) cos(x+) sin(x-) (14) 



can also be non-zero. We discuss the meaning of this new order parameter, which is non- 
local in the growth-direction spatial coordinate, at greater length later. Finite expectation 
values for either operator require broken symmetry since the Hamiltonian is invariant under 
rotations about the x axis in spin space. 

Given the density matrix (|I2D, the Hartree-Fock single particle Hamiltonian is readily 
evaluated: 

2h HF = -F s + (AH -F s - A v ) cos(#) cos( X+ ) cos( X -)r 2 - (F s + A z ) sin(0) sm( X -)a x 

- (F D + A t ) cos(#) sin(x+) cos(x-)t x + F D cos(#) cos(x+) sm( X -)(r z T x 

- F D sin(6)(cos(x+) - cos(x-))a x T x - F s cos(6)sm(x + )sm(x-)(T z T z , (15) 

where F$ = F + + F_ and Fjj = F + — F— When the variational parameters are at an ex- 
tremum of the Hartree-Fock energy functional (|TTD, z 1 and z 2 are eigenvectors of matrix flT5|) 
and if the extremum is the minimum these eigenvectors will correspond to the two lowest en- 
ergy eigenvalues, E± F and E% F . The Hartree-Fock approximation for the charged excitation 
energy gap, measured by quantum Hall transport experiments, is E gap = E^ F — E^ F . 

III. THE VANISHING ZEEMAN COUPLING LIMIT 

At A z = 0, we expect the ground state spin-magnetization directions in the two quan- 
tum wells to be in opposition. The collinear spin configuration corresponds to taking the 
parameter 6 to be zero in the ansatz (|10|). This reduces the number of free parameters in the 
variational wavefunction to two, namely, the polar angles for the up-spin and down-spin pseu- 
dospin orientations. With 8 = 0, the spinors (fL0| ) reduce to z 1 = (cos(x^/2), 0, sin(x^/2), 0) 
and z 2 = (0,cos(xV 2 ),0,sin(xV 2 )) ; and the 

energy per orbital state (|TT| ) reads 

e = -A y cos(x+) cos(x_) + 2H cos 2 (x+) cos 2 (x-) 
- A t sin(x+) cos(x_) - F+ - F_(cos 2 (x+) cos 2 (x_) + sin 2 (x+) sin 2 (x-)) . (16) 

The first two terms on the right hand side of this equation represent the electrostatic 
energy which is minimized when the 'flat band' condition, cos(x + ) cos(x~) = Av /AH, is 
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satisfied. If these were the only terms in the energy expression, the charge transfer between 
layers would perfectly screen the external interlayer electric field parameterized here by Ay. 
Since we have not included Zeeman coupling, any difference between the pseudospin angles 
of up-spins and down-spins, i.e. X- 7^ 0, represents a broken symmetry. The exchange 
energy is maximized in magnitude by localizing both up-spins and down-spins in one of the 
layers; at Ay = this objective can be made consistent with minimizing the electrostatic 
energy by choosing opposite layers for the two spins, i.e., x+ = 7r /2- Minimizing with 
respect to X- then gives cos(x-) = A t /2F_ for A t < 2F_. Otherwise cos(x-) = 1 and 
there is no broken symmetry. The spin- up and spin-down pseudospins for A z = 0, Ay = 
are schematically illustrated in Fig. |2]. Here and in all figures below, we plot results for 
layer separation d = 1 in units of the magnetic length i yielding the following interaction 
strengths: H = 0.5, F + = 0.9545, and F_ = 0.2988, where energy is measured in units of 
e 2 /e£. The mean pseudospin orientation is in the x direction which corresponds to equal 
distribution of charge between the layers. For A t > 2F_ the two pseudospins are parallel 
(X- — 0) an d the spin polarization in each well is zero, as also indicated in Fig. ^. For 
A t < 2F_, z 1 and z 2 have pseudospin orientations tilted from the x direction upwards and 
downwards, respectively, x~ 7^ and, hence, the order parameter O zz is non-zero. As the tilt 
angle increases with decreasing A t , opposite spin polarizations builds up in the two layers 
reaching a maximum as A t — ► 0. Then the two pseudospins are in opposition and aligned in 
the z direction. The non-zero value of O zz in the ordered region reflects the broken SU(2) 
symmetry (full rotational symmetry) of spin space. Note that the order parameter O xz is 
zero along the A z = 0, Ay = line since x+ — 7r /2- 

We now turn our attention to A z = 0, At = line in the phase diagram. In this case the 
energy is minimized when one of the spinors, e.g. z 1 , is localized in the top well, i.e., %t = 
so that X- = X+ = Xi- Minimizing the energy functional ([16]) with respect to X- gives 

for 2F_ < Ay < AH - 2F_, cos 2 (x_) = for Ay < 2F_ and cos 2 (x_) = 1 for Ay > AH - 
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2F_. For bias potential Ay < 2F_ the up and down spin orbitals are localized in opposite 
wells and no charge transfer is produced, as illustrated in Fig. |2|. For Ay > AH — 2F__, the 
up and down spin orbitals are localized in the top well and the charge transfer is complete. 
At intermediate values of Ay, Eq. ([T7|) applies. In the interval 2F_ < Ay < AH — 2F_, both 
O zz 7^ and O xz ^ 0. The non-zero value of O zz reflects broken SU(2) spin-symmetry in the 
ground state while the non-zero value of O xz at A t = reflects a U(l) broken pseudospin 
symmetry. At At = the Hamiltonian is invariant under rotations about the z axis in 
pseudospin space; the number of particles in each layer is a good quantum number. This 
symmetry is broken since the orbital states with variational spinor z 2 possess interlayer 
phase coherence (see Fig. 0). As in the v = 1 case, spontaneous interlayer phase coherence 
will lead to a variety of interesting effects^! which we do not pursue further here. We return 
to this limit again in Section V. 

We have been unable to obtain general analytic expression for the minimum energy state 
for both Ay and A t non-zero. Instead we define C = cos(x+) cos(x_), S = sin(% + ) sin(x_) 
and expand around the disordered state for which 5 = 0. We consider S as an order 
parameter and perform a Ginzburg-Landau analysis. 

e = e (C) + t 2 {C)S 2 + \t A {C)S A + .... (18) 

Comparing with Eq. (0), we obtain the following results for the coefficients of this Ginzburg 
Landau expansion: 



e (C) = -AyC + 2HC 2 - A ty J(l - C 2 ) - F+ - F^C 2 



e 2 (C) 



2(1 - C 2 fl 2 



A 1 _|_ AC 2 

^ = TJT^cW (19) 

In the normal state, S is zero and C is determined my minimizing eo(C), implying the 
following relationship between Ay and C: 

Ay = (AH - 2F_)C + A t (i _^ 2)1/2 (20) 
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The normal state becomes unstable when ^{C) becomes negative; on the phase boundary, 
C equals 

C* = (1 - (A t /2F_) 2 / 3 ) 1/2 - (21) 

As the bias potential Ay increases, 62(C) increases, crossing zero at the critical bias potential 

A* v = [AH - 2F_ + A 2/3 (2F_) 1/3 ](1 - (At/2F_) 2/3 ) 1/2 . (22) 

Eq. (|22"D describes the projection of the mean-field phase boundary onto the A 2 = plane, 
illustrated in Fig. Q. 

Among the observables which are singular at this T = phase transition, the most 
experimentally accessi ble§ is the internal differential capacitance of the bilayer system: 

In the normal state, an equation for Ci nt can be obtained by differentiating both sides of 
Eq. ( p0|) with respect to Ay: 

9 

int ~ AH-2F_ + A t (l -C 2 )^/ 2 1 1 

Only the first term in the denominator appears in the electrostatic approximation for this 
capacitance. A similar expression, valid near the phase boundary, may be derived for the 
ordered state by minimizing the Ginzburg-Landau energy expression with respect to S. We 
find that 

e w e (C) - A(C* - C) 2 . (25) 

where 

_ 18F_(l-(A f /2F_) 2 / 3 ) 
A " 5-4(A,/2F_)2/3 ' ^ 

The additional condensation contribution to the energy alters the inverse capacitance. Dif- 
ferentiating Eq. ( P5D first with respect to C and then with respect to Ay we find that as 
the phase boundary is approached from the ordered state side 
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« order 
int 



2 



(27) 



AH - 2F_ + A t (l - C 2 )- 3 / 2 - A ' 



In mean-field-theory the capacitance has a jump discontinuity at the phase boundary. This 
observable is the analog, for this bias tuned quantum phase transition, of the specific heat at 
a temperature tuned phase transition. The differential capacitance is singular at all points 
in the two-dimensional boundary of the ordered state region in the three-dimensional phase 
space. This singular behavior should be accessible to experiment and may be interesting 
to study, especially since the nature of the broken symmetry in the ordered state crosses 
over from SU(2) to U(l) on going away from the A z = plane. The mean-file critical 
exponent will be replaced by the specific heat exponent of a three-dimensional classical 
system, presumably that of the 3D Heisenberg model at A z = and that of the 3D XY 
model at A z ^ 0. In Section V we give an analytical expression for the capacitance jump 
for finite Zeeman coupling in the limit of vanishing interlayer tunneling. 



The Hartree-Fock equations for the zero bias potential limit have been solved previously 
by Das Sarma et aW in their pioneering work on v = 2 broken symmetry ground states. In 
this section we report analytic solutions of the Hartree-Fock equations which provide some 
additional insight. When Ay = 0, we expect the energy minimum to occur for equal charge 
density in the two layers, i.e., at (t z ) = 2 cos(6 l ) cos(x+) cos(x_) = 0. It is clear that the 
energy is minimized by x+ — 7r/2, since this choice does not frustrate any of the other energy 
terms. The Hartree energy then vanishes and the total energy expression fllTD simplifies to 



IV. THE VANISHING BIAS POTENTIAL LIMIT 



e = — 



At cos(0) cos(x-) — A z sin(0) sin(x-) — F + — F_ sin 2 (%_). 



(28) 



Minimizing with respect to 6 we find that at the extremum 




(29) 



Minimizing with respect to X- an d using Eq. (|2~9"|) we find that 
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(A f - A 2 /A t ) cos(fl) = 2F_ cos(x-). (30) 
These two equations can be combined to find explicit equations for cos(#) and sin(x_): 

cos 2 (£) = ^K A? " A ' )2 " ^ 2F - A ^ 



(A? - Al) 



,3 



2 , N 4A 2 F 2 — (A? — A 2 ) 2 , . 

Broken symmetry states occur when both of these equations have solutions. For A t 
larger than 2F_ and A z smaller than 



A t (A t - 2F_) (32) 

the expression for sin 2 (x_) is negative. The minimum energy then occurs for sin(x_) = 
and 6 = and the variational spinors reduce to the up-spin and down-spin symmetric 
states, as illustrated in Fig. [| These states are expected to have the lowest two energies for 
non- interacting electrons for A z < A t . For A z larger than 



F 2 + A 2 - F_ (33) 

the expression for cos 2 (#) is negative. The minimum energy then occurs for 9 = n/2, xt = 0, 
and xi = ft- I n this case the variational spinors are both positive eigenvalue eigenstates of 
a x . The wavefunction in this region is the non- interacting electron state for A z > A t . In 
the ordered region, O zz is non-zero and the spin U(l) symmetry is broken. The magnetic 
configuration (see Fig. [|) is that of a canted antiferromagnet as observed by Zheng et a/..i 
Each of the two spinors z 1 and z 2 has non-zero spin up and spin down components in the 
ordered region for A z > 0. However, the weights of spin down component of z 1 and spin 
up component of z 2 decrease with decreasing Zeeman coupling and the state continuously 
develops into the one illustrated in Fig. [| in the limit of A z = 0. Note that A™ ax cx A 2 
at small A t so the broken symmetry state is not robust against Zeeman coupling when 
A t is small. Note also that when A t is large A™ m = A t — F_ — F 2 /2A t + . . . while 
A™ ax = A t — .F_ + F 2 /2A t + . . .; the broken symmetry state interpolates between the fully 
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spin-polarized and fully pseudospin polarized states and occurs only over a narrow interval 
of Zeeman fields surrounding A z = A t — F_. The field at which the crossover occurs is 
shifted from A t because of the exchange energy is more favorable for the full spin-polarized 
state than for the fully pseudospin-polarized state. 



V. THE VANISHING INTERLAYER TUNNELING LIMIT 

For A t = 0, we expect x+ — 0. Numerical calculations presented in the next section verify 
that this choice of x+ gives the optimal variational energy for vanishing interlayer tunneling. 
Given x+ — 0, it can be seen that the energy is minimized by 9 = X-i leaving a single 
parameter to be determined by minimizing — Ayx — A 2 (l — x)— F + — — x) 2 + x 2 ) + 2Hx 2 
with respect to x — cos 2 (9) = cos 2 (x-)- We find that cos(#) = cos(x-) = for Ay < 
2F_ + A z , cos(9) = cos(x-) = 1 for A v > AH - 2F_ + A z , and that 

2(n\ Ay - A z - 2F_ 

cos (6) = a{h _ f ) (34) 

in the intervening interval. The low bias state is fully spin polarized along the Zeeman 
field direction and has the charge equally distributed between the two layers, as indicated 
in Fig. |j. The high bias state is spin unpolarized and has all electrons in the top layer. For 
non-interacting electrons, the transition between these two states is first order and occurs at 
Ay = A z . In the Hartree approximation, charge transfer occurs linearly with bias voltage, 
since the double layer system acts like a capacitor. In this approximation, the minority spin 
Landau level in the top layer and the majority spin Landau level in the bottom layer are 
degenerate throughout the interval of bias voltages where charge transfer takes place. We 
show below that an energy gap exists throughout this interval, and that the gap is supported 
by the development of a broken symmetry state. However, the broken symmetry is not the 
canted antiferromagnet order of the Ay = plane but rather something closely akin to the 
spontaneous interlayer phase coherence which occurs at v — 1 in double-layer systems. The 
projection of the phase diagram onto the A t = plane is illustrated in Fig. f|. 
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In the ordered portion of the A t = plane, O zz = and 

°- = > ~ ( %% 2H ^ 12 - < 35 > 

In the ordered state, the Landau level orbitals with spinors z l and z 2 spontaneously develop 
interlayer phase coherence and spin-polarization in the z direction. (The Zeeman field is ori- 
ented in the x direction.) One Landau level orbital has up and down spin projections whose 
pseudospin orientations are proportional to [sin(%_), 0, cos(x-)] and [sin(x_), 0, — cos(x_)] 
respectively; the spin projections of the other Landau level orbital have opposing pseudospin 
orientation and altered weights (see Fig. ^). Spontaneous interlayer phase coherence corre- 
sponds to non-zero pseudospin components in the x direction. The occurrence of this type 
of order is driven by interlayer interactions; the spatial correlations between electrons in 
different layers are improved when it is established. 

We note that the solutions we have found for the A t — > limit of the A z = plane and 
the A z — ► limit of the A t = plane differ. Evidently, the order of limits matters in this 
instance. It turns out that when both these terms are absent, energy minimization fixes 
a value only for the product cos(#) cos(x + ) and not for these factors individually. In the 
A t = plane, the state is invariant under a simultaneous rotation of all pseudospins around 
the z axis and under a rotations of spins around the x axis. These invariances are explicitly 
accounted for in limiting the number of parameters in our variational wavefunction. For 
A z = 0, the state is invariant under rotation about any spin-axis. When both external 
fields are set to zero, all symmetries are present, leading to additional soft modes and a 
set of variational wavefunctions which have the same energy. When these field parameters 
approach zero, O zz ^> O xz if A t ^> A z , and O xz ^> O zz if A z ^> A t . From the above solution 
for the variation parameters, it is possible to show that for A t = 0, E gap = Fp, a constant 
throughout the ordered region of the phase diagram. The gap arises entirely because of the 
interlayer interactions which give rise to spontaneous interlayer phase coherence. 

In the A t = plane, (t z ) = 2cos 2 (6 l ). Using Eq. (^) we can derive an analytical 
expression for the internal capacitance in the ordered phase: 
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In the normal state, no charge transfer is produced by changing Ay in either low or high 

bias regions of the phase diagram. The internal capacitance is hence zero throughout the 

normal state region and has a jump, given by (^), along the Ay = A z + 2F_ and Ay = 

A z + AH - 2F_ lines. 



VI. NUMERICAL RESULTS 

We supplement the analytical results discussed in the previous sections with results 
obtained from numerical solutions of the four-dimensional Hartree-Fock equation outlined 
in Section II. An overview on the order-disorder physics of bilayer quantum Hall systems 
at v = 2 is provided in Fig. [5|, a three-dimensional phase diagram whose axes are the three 
one-body external fields, Ay, At, and A z . In this diagram, the region with spontaneous 
order is enclosed by two surfaces. In the limit of no interactions these collapse onto the 
single surface, Ay + A^ = A 2 Z , where the 2nd and 3rd Landau levels are degenerate. Broken 
symmetry states can occur at integer filling factors when interactions dominate the energetic 
splitting of nearby Landau levels. The width of the ordered region tends to expand when 
the single-particle field parameters are small and all four single-particle Landau levels are 
close to degeneracy. The intersection of this phase diagram with the Ay = plane recovers 
the results obtained originally by Das Sarma et al. and derived analytically in Section IV. 
In the ordered region within the Ay = plane, sin(x + ) = 1, and O zz ^ while O xz = 0. 
When A t ^ 0, inter layer phase coherence is not spontaneous so O zz ^ represents only the 
£-axis spin-rotation (U(l)) broken symmetry. In the ordered portion of the A t = plane 
and for A z > 0, O xz ^ 0, but O zz = 0. The order here involves spontaneous interlayer phase 
coherence and spin order; as at v = 1 we can expect interesting physics to occur for small A t 
in tilted fields in this region of the phase diagram. The ordered state has x-axis spin-rotation 
and z-axis pseudospin-rotation broken symmetries ( U(1)®U(1) broken symmetry) so we 
expect two distinct Goldstone modes to appear. At A z = the broken symmetry changes 
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to SU(2)(g)U(l) in the ordered region of A t = plane. The phase diagram intersections 
with the A z = and A t = planes are in agreement with the analytic results presented in 
Sections III and V respectively. In the ordered region of the A 2 = plane, both O zz and O xz 
are non-zero and are indices of the same broken symmetry. In this case, the Hamiltonian is 
spin-rotationally invariant so we have a SU(2) broken symmetry. 

In experimental samples, Ay and At can easily be made large but A z is usually small. 
Since it is easy to tune the bias voltage Ay over large ranges experimentally, we examine be- 
low the dependence of physical observables on this parameter with the other two parameters 
held fixed. 

In Figs. 0(a)- (c) we plot the bias dependence of the following quantities, all of which 
are experimentally accessible: the two order parameters O zz and O xz , the Hartree-Fock 
energy gap E gap , the difference between top layer and bottom layer charge densities (r 2 ), 
and the differential capacitance of the bilayer system, d(r z ) /dAy. The bias dependence 
of the energy gap has been extensively stud.ed experimentally^ we expect the present 
Hartree-Fock theory to capture the main trends^ in these observations although we know 
that spin-texture excitations can be important^ and that disorder, not accounted for here, 
plays an important role in experimental samples. The most interesting quantity in our view 
is the internal differential capacitance, which shows singular behavior at the phase boundary. 
Results are plotted here for (A t , A z ) = (0,0.02), (0.02,0.02), and (0.1,0.02). These choices 
are motivated by the fact that it is difficult to fabricate samples with values of the Zeeman 
coupling in interaction units e 2 /e£ which differ from 0.02 by more than a factor of two at 
v = 2, while it is relatively easy to vary A t from a large value to one which is immeasurably 
small, by changing the Aluminum content in the barrier separating the quantum wells. 

Numerical calculations for A t = and A z = 0.02, presented in Fig. 0(a), confirm 
analytical results of Section V: the order parameter O zz = while O xz is non-zero for 
A v > A z + 2F_ = 0.6176 and Ay < A z + AH - 2F_ = 1.4024; the quasiparticle excitation 
gap E gap is constant in the ordered region; the charge is equally distributed between layers 
in the low Ay disordered region, all electrons are in the top layer in the high Ay disordered 
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region and the charge is transferred between wells linearly with Ay in the ordered region. In 
Fig. ^](b) we show numerical results for A t = A z = 0.02. Both low and high Ay disordered 
regions are still present for this choice of external fields. The ordered region with both O xz 
and O zz non-zero is slightly shifted to lower bias potentials and the step in Ci nt is reduced at 
the right phase boundary. Fig. ^|(c) corresponds to large enough tunneling energy (A t = 0.1) 
so that the Ay = point falls into the ordered region of the zero-bias plane of the phase 
diagram (see Fig. [3]). Hence, only the high Ay phase boundary is present and is shifted 
further to lower Ay, qualitatively consistent with our previous analysis of the A z = phase 
plane. 

VII. DISCUSSION 

The variational ansatz (|K]) on which we have based the analytic portion our analysis 
carries only 3 parameters. One might wonder if more general wavefunctions with more 
parameters could unearth lower energy states. Let us consider how many parameters at 
most can be involved in the most general translationally invariant LLL state of bilayer 
electrons with spin, at a filling v = 2. Such a state will correspond, at each spatial point or 
LLL orbital, to some 2-fermion state where each fermion is described by a 4-spinor. More 
generally, consider N fermions, each a d-spinor. There are d\/(N\(d — N)\) combinations 
for placing N fermions in d states. Each such combination corresponds to an independent 
N-fermion state. These form a basis and a general state will be a linear combination of these 
with complex coefficients . Since each complex coefficient carries two real parameters the 
total number of parameters in the most general state will be n = 2d\/(N\(d — N)\) — 2, 
where 2 parameters have been subtracted out which respectively account for normalisation 
and an overall phase. For our problem where d=4 and N=2, n = 10; much larger than the 
number (three) of parameters in our trial spinors. 

However not all these n independent states can be written as a single outer product of 
two 4-spinors, as is required for the Slater determinant form (|IJ) of the Hartree-Fock trial 
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states. For example , the linear combination of 2-electron states 



\TITI) + \BtBl) 



(37) 



(where T,B refer to the two layers) cannot be written as an outer product of two 4-spinors 
of the form 



( a x \ 
\a 4 J 



h 



(38) 



for any values of the constants and hi which would need to satisfy the conflicting conditions 
ai&2 = ^3^4 = 1; while = 0362 = 0. Whether such variational states which, for each 
spatial point, are sums of products of spinors (as for instance in fl37|)) can be candidates for 
the ground state of different phases is an interesting question. Apriori they are not ruled 
out by any basic superselection rules. However, broken symmetry states with which we are 
familiar, for example superconducting and magnetic states, can be described using an order 
parameter field functional integral language in which a saddle point approximation to the 
order parameter functional integral corresponds to a Hartree-Fock approximation for the 
broken symmetry ground state. We assume that to be the case for double-layer quantum 
Hall states as welfl. In particular, we believe that the set of possible broken symmetries 
can be identified by looking only at single Slater determinant states. The positions of 
the phase boundaries will be altered by quantum and, at finite temperature, also thermal 
fluctuations in the order parameter field. We therefore are motivated to identify the number 
of parameters needed to characterize all possible Hartree-Fock states. 

A product of N orthonormal d-plets will carry Nd complex parameters to start with. 
This will be reduced by N normalization constraints, N(N — l)/2 complex orthogonality 
constraints, and iV 2 symmetry constraints since unitary transformations among the N the 
single-particle states will not yield a distinct N-fermion state. This leaves behind n = 
2dN - N-2N(N-l)/2 - N 2 = 2N(d - N) real parameters. Notice that this is just the 
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dimensionality of the coset space U (d) / (U(N) x U(d — N)). (See Arovas et a]0 who also do 
this counting). For d — 4 and N = 2 the number of free parameters in Hartree-Fock states 
is n = 8; for N — 1, the case relevant to double-layer systems at v — 1, n — 6. For our 
problem, the number of free parameters is further reduced by symmetries of the interaction 
Hamiltonian which is invariant under arbitrary (3-parameter) rotations in spin-space and 
under (1 more parameter) rotations about the z axis in pseudospin space. Of course these 
symmetries are partially broken by one-body terms in the Hamiltonian. However, such 
breaking simply selects a known preferred direction for the total spin-polarization (along 
the Zeeman coupling direction ) and for the projection of the total pseudospin polarization 
onto the x — y plane (the x direction.) The end result is that the number of free parameters 
is 4 for v = 2. We have used 3 parameters based on the assumption that in the grounds 
state pseudospin projections in each layer will be in the x — z plane. 



VIII. CONCLUDING REMARK 

The richness of order-disorder physics in double-layer quantum Hall systems at v = 2 
suggests that this should be an inviting system for experiment. To date, only transport!^ 
and optical studies, which probe the state rather indirectly have been completed. We be- 
lieve that capacitance studies of double-layer systems as a function of an interlayer bias 
potential will prove more interesting. In particular we expect the differential capacitance, 
an experimentally accessible quantity which is the analog for these zero-temperature phase 
transitions of the specific heat in thermal phase transitions, to be a very useful probe in 
studying changes in singular dependence on bias potential with temperature, Zeeman cou- 
pling strength, and interlayer tunneling amplitude. 
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FIGURES 

FIG. 1. Schematic of the conduction band edge profile and one-particle energy levels in a 
biased double quantum well structure. The four one-particle states are products of orbitals split 
due to both interlayer tunneling and a bias potential, and spinors split by the Zeeman coupling. 
Interactions become important and broken symmetry states occur at v = 2 when the second and 
third levels, and sometimes also the first and fourth levels, are close to degeneracy 

FIG. 2. A z = plane of the phase diagram. The ground state has SU(2) broken symmetry 
throughout the shaded region and U(1)®SU(2) broken symmetry along the darkened portion of the 
A z = — At = line. The dashed-line arrows labeled (i, a) illustrate schematically the pseudospin 
of spin-cr projections of the spinor z l at a given point in the phase diagram. An upward arrow 
represents an orbital localized in the top layer, downward arrow an orbital in the bottom layer, while 
a finite tilt from the vertical direction indicates interlayer phase coherence. The solid-line arrows 
show schematically the orientation of the spin polarization in top (upper box) and bottom (lower 
box) quantum wells. The length of the arrow represents the magnitude of the spin polarization in 
that layer while its thickness represents the fraction of the charge density. 

FIG. 3. Ay = plane of the phase diagram. The ground state has a U(l) broken spin sym- 
metry in the shaded region. The dashed-lines and solid-lines summarize respectively the spin and 
pseudospin structures of the ground state as in Fig. 2. 

FIG. 4. At = plane of the phase diagram. The dashed-lines and solid-lines summarize 
respectively the spin and pseudospin structures of the ground state as in Fig. 2. Spontaneous 
interlayer phase coherence is developed in the shaded region. 

FIG. 5. 3D phase diagram for bilayer system at v = 2. Ordered states occur in the volume 
enclosed by the shaded and transparent surfaces. 
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FIG. 6. Order parameters, quasiparticle excitation gap, charge distribution, and internal ca- 
pacitance as a function of bias potential for A t = and A z = 0.02 (a), A t = A z = 0.02 (b), and 
A t = 0.1 and A z = 0.02 (c). 
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